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^ ' Abstract. With a mean field approach, the heavy Fermi liquid in the two-dimensional 

C/3 , Kondo lattice model is carefully considered in the presence of short-range antiferromagnetic 

correlations. As the ratio of the local Heisenberg superexchange coupling to the Kondo coupling 
increases, the Fermi surface structure changes dramatically. From the analysis of the ground 
state energy density, multiple Lifshitz type phase transitions occur at zero temperature. 
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1. Introduction 

Quantum phase transitions have attracted much interest in studying correlated electron systems. 
^ ■ An electronic phase transition associated with the change of Fermi surface (FS) topology, the 

•/^ ', so-called Lifshitz transitional], can be induced without any spontaneous symmetry breaking and 

' local order parameter. The Kondo lattice model is a prototype model and believed to capture 

. the basic physics of heavy fermion materials. The huge mass enhancement of the quasiparticles 

can be attributed to the coherent superposition of individual Kondo screening clouds, and the 
resulting metallic state is characterized by a large FS with the Luttinger volume containing both 
conduction electrons and localized moments. Competing with the Kondo singlet formation, the 
localized spins indirectly interact with each other via magnetic polarization of the conduction 
electrons - the Ruderman-Kittel-Kasuya-Yosida interaction. Such an interaction dominates at 
^ I low values of the Kondo exchange coupling and is the driving force for the antiferromagnetic 

^ ■ (AFM) long-range order and quantum phase transitions [21 [3]. So far most of investigations focus 

on the possible FS reconstruction around the magnetic quantum critical point [H El El [3 [8]. 
However, the FS topology in the paramagnetic heavy fermi liquid phase may also be drastically 
changed by the short-range AFM spin correlations between the localized spins, leading to the 
Lifshitz phase transitions|9]. The nature of such a quantum phase transition has not been 
thoroughly explored yet. 

In this paper, we consider the two-dimensional Kondo lattice model with the Heisenberg 
AFM superexchange coupling between localized spins. By introducing uniform short-range 
AFM valence-bond and Kondo screening parameters, a fermionic mean-field theory is derived 
and carefully re-examined. Away from half-filling, at the conduction electron density ric = 0.85, 
for example, the possible changes of FS topology in the paramagnetic heavy Fermi liquid phase 
are considered carefully as increasing the short-range AFM spin correlations. 



2. Mean field theory 

The model Hamiltonian defined on a square lattice is given by 

H = ekCkfjCkCT + JR^Si ■ Si + Jn^Si ■ Sj. (1) 

k,cr i (ij) 

The spin-1/2 operators of the local magnetic moments have the fermionic representation Sj = 

2 X^o-cr' fia'^cTCT' fia' with a local constraint X^o- fiafio^ — where r is the Pauli matrices. Following 
the large- A'^ fermionic approach [TU], the Kondo spin exchange and Heisenberg superexchange 
terms can be expressed up to a chemical potential shift as 

Sj • Sj = — - ^ ficrfjafja-'fia', Sj • Sj = — - ^ firjCjaC^j^'fia'-, (2) 

<T(t' aa' 

then uniform short-range AFM valence bond and Kondo screening order parameters can be 
introduced as 

x = -E(/U+'-)' ^ = E(4/-)- (3) 



To avoid the incidental degeneracy of the conduction electron band on a square lattice, we 
choose ek = —2t (cos + cos ky) + At' cos kx cos ky — fi, where t and t' are the first and second 
nearest neighbor hoping matrix elements, respectively, while /i is the chemical potential, which 
should be determined self-consistently by the density of the conduction electrons ric- Under the 
uniform mean-field approximation, the f-spinons form a very narrow band with the dispersion 
Xk = JhX (cos kx + cos ky) + A, where A is the Lagrangian multiplier to be used to impose the 
local constraint on average. 

Thus the corresponding mean-field Hamiltonian reads 

with Eq = N {—X + JhX^ + JrV'^ 1'^) ■ The quasiparticle excitation spectra can be easily 
obtained 

^k =^[(ek + Xk)±H^k], (5) 

which implies that the conduction electron band ek has a finite hybridization with the spinon 

band Xk- Here W\f_ = \J (ek — Xk)^ + {JRy)^ ■ Accordingly, the ground-state energy density can 
be evaluated as 

= I E ^k ^ (-^k ) - A + JhX^ + \JrV\ (6) 

k,± 

where 9{—E\^) is the theta function. Then the self-consistent equations for the mean-field 
variables and A and the chemical potential fi can be deduced from the relations 

^-0 ^-0 ^-0 n --^ (7) 
3. Results and Discussions 

In the following, we will assume that t'/t = 0.1 and ric = 0.85. When the self-consistent 
equations are carefully evaluated, we find that the mean-field AFM order parameter x is always 
positive in the range < Jh/Jr < 3 so that the resulting state is a stable paramagnetic metal. 



In Fig.l, as the strength of the AFM spin fluctuations grows up, we present the evolution 
of the band structure of the renormalized heavy quasiparticles around the Fermi level in the 
direction (0,0) (7r,vr) — > (vr, 0) — > (0,0) of the first Brillouin zone. Fig. la to Fig.lh 
correspond to Jx/t = 2.0 and x = Jh/Jk = 0, 0.11, 0.12, 0.234, 0.235, 1.2, 2.76, and 2.77, 
respectively. At the critical values x = 0.1181, the point (vr, vr) changes from the local maximum 
to local minimum of the quasiparticle band, while at a; = 0.2341, the separated band around 
the momentum (vr, 7r/2) starts to move away from the Fermi level. As the ratio of the coupling 
strengths increases, the separated band moves to the momentum (7r, 0), and reaches the Fermi 
level again at x = 2.77. 



(a) x=0 



(b)x=0.11 



-0.050 




0.025 
0.000 

'uj 

-0.025 
-0.050 
0.025 
0.000 
■'"'.0.025 



(e) x=0.235 



(g) x=2.76 



-0.0501 

(0,0) 



(f)x=1.20 



(h) x=2.77 



(1,0) 



(0,0) (0,0) 



{n,n) 



(0,0) 



Figure 1. The lower renormalized quasipar- 
ticle band in the direction (0, 0) — )• (vr, vr) — >■ 
(vr, 0) — )• (0, 0) as increasing the strength of 
the AFM spin fluctuations. 
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Figure 2. The first-order derivative of the 
ground state energy as a function of the 
parameter x = Jh/Jk- There are three 
singular points, Xc,i = 0.1181, Xc,2 = 0.2341, 
and Xc,3 = 2.77. The first one corresponds 
to a first order transition, while the latter 
two correspond to two second order phase 
transitions, respectively. 



From the quantum phase transition aspects, we calculate the ground state energy density 
Eg and its first-order derivative with respect to the ratio of the coupling parameters x. The 
numerical results are displayed in Fig. 2. We find that there are three non-analytical points. £g 
is finite and continuous in the parameter range < x < 3, However, its first-order derivative 
has a large jump at xic = 0.1181, corresponding to a first-order quantum phase transition. 
Moreover, two small kinks appear at X2c = 0.2341 and xsc = 2.77 in the first-order derivative, 
which correspond to the jumps in the second-order derivative of Eg. So X2c and xsc denote two 
second-order quantum phase transitions. 

Once the renormalized quasiparticle band structure is available, the corresponding FS can 
be easily obtained. The obtained FS is the large hole-like one and shown as the shaded area. 
In Fig.3a-3g, the center of the FS is shifted from (0,0) to (vr,vr). For a fixed Jx/t = 2.0, we 
can see that the FS is a hole-like circle around (vr,vr) for the parameter range < x < 0.11. 
At xic = 0.1181, the topology of the FS starts to change: a small circle emerges in the center 
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Figure 3. The structure of the FS changes as increasing the strength of the AFM spin 
fluctuations. The shaded area represents the hole-hke FS. The center of the FS is shifted at 
(7r,7r) in (a)-(g). 



of the deformed large square FS. As x is further increased, both circles expand and the small 
one is deformed into a rotated square. Up to X2c = 0.2341, the two deformed circles intersect 
each other and then decompose into four Fermi pockets. When x increases to x^c = 2.77, the 
FS pockets are reconnected again and form two closed square-like circles. Three critical values 
correspond to three different quantum phase transitions, which belong to the category of Lifshitz 
phase transitions. 

To some extent our present mean-field theory captures the heavy-fermion liquid physics of 
the Kondo-Heisenberg lattice systems, especially the Fermi surface evolution of the renormalized 
heavy quasiparticles as the short-range AFM spin correlations between the localized magnetic 
moments are gradually increased. In order to put the present results on a more solid ground, 
further investigation beyond the mean-field theory is certainly needed. 
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